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We study classically unstable string type configurations and compute the renormalized vacuum 
polarization energies that arise from fermion fluctuations in a 2+1 dimensional analog of the stan- 
dard model. We then search for a minimum of the total energy (classical plus vacuum polarization 
energies) by varying the profile functions that characterize the string. We find that typical string con- 
figurations bind numerous fermions and that populating these levels is beneficial to further decrease 
the total energy. Ultimately our goal is to explore the stabilization of string type configurations in 
the standard model through quantum effects. 

We compute the vacuum polarization energy within the phase shift formalism which identifies 
terms in the Born series for scattering data and Feynman diagrams. This approach allows us to 
implement standard renormalization conditions of perturbation theory and thus yields the unam- 
biguous result for this non-perturbative contribution to the total energy. 

PACS numbers: 03.65Sq, 03.70+k, 11.27+d 



I. INTRODUCTION 



In the past decades the perturbative treatment of the electroweak standard model has proven to be a very powerful 
tool to describe the properties and interactions of elementary particles within a wide energy regime. On the other 
hand the role and even the existence of non-perturbative solutions in this model is still quite uncertain. While the 
standard model does not contain topological solitons, one can construct nontopological string solutions to the classical 
equations of motion, called Z-strings or electroweak strings Q, |^ |^ . In the absence of topological arguments, however, 
one is not guaranteed that these classical solutions actually correspond to true local minima of the full effective energy, 
or even the classical energy. Indeed, Naculich Q has shown that in the limit of weak coupling, fermion fluctuations 
destabilize the string solution. This analysis leaves open the possibility, however, that deformed string solutions 
could exist, which would be local minima of the full effective energy. Since fermions are tightly bound in the string 
background, one potential mechanism for restoring stability is for fermions to bind to the string, yielding a lower total 
energy than a corresponding density of free fermions. When including this effect, however, one must also take into 
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account the shifts in the zero-point energies of all the unoccupied modes as well, computed consistently in a standard 
renormalization scheme. 

If stabilized by the fermion binding mechanism, electroweak strings could have significant cosmological consequences 
,2,] . A network of strings could contribute to the dark energy that is required to explain the recentl y ob served cosmic 
acceleration [Tj, Ig- However, a complete dynamical description of that scenario is still missing ahOlhll. Also, as 
pointed out by Nambu 0, Z-strings are expected to terminate in monopole-antimonopole pairs, which could give 
rise to a primordial magnetic field. Furthermore, a network of stable strings at the electroweak phase transition 
would provide a scenario for electroweak baryogenesis without requiring a first-order phase transition |12|. The 
strings provide out-of-equilibrium regions, and the core of the string has copious baryon number violation due to the 
suppressed Higgs condensate. They thus provide an alternative to the usual idea of bubble-nucleation baryogenesis, 
which requires a first order phase transition to go out of thermal equilibrium. 

From a theoretical point of view, the quan tum properties of Z-strings have been connected to non-perturbative 
anomalies |0. Decoupling arguments [T^llSf suggest that when a fermion's mass is made very large by increasing its 
Yukawa coupling, soliton configurations should appear in the low-energy spectrum to maintain cancellation of such 
anomalies. In these calculations, the analysis of the fermion determinant - whose logarithm is the sum over zero-point 
energies - is essential. 

A first attempt at a full calculation of the quantum corrections to the Z-string energy was carried out in 
Those authors were only able to compare the energies of two string configurations, rather than comparing a single 
string configuration to the vacuum. Furthermore, they used a proper time prescription, in which it was not possible 
to explicitly split off the divergent parts from the fermion determinant. Thus the final answer was expressed at large 
but finite cut-off in terms of two large terms that could only be computed numerically, representing the fermion 
determinant and the counterterms. In this formulation, no numerically stable results could be extracted for the limit 
"cut-off to infinity." In the present work, wc will employ phase shift techniques |17| . which allow us to make this 
separation cleanly and unambiguously. 

In this paper we restrict ourselves to the case of 2 -I- 1 dimensions. We make this restriction to simplify the 
renormalization procedure. However, from the experience gathered in the case of QED magnetic flux tubes it is 
reasonable to exp ect that the qualitative behavior might be very similar to the actual 3-1-1 dimensional theory. The 
techniques of allow for a straightforward generalization of the calculational procedure used here to that case as well. 
We also consider only the weak interactions, neglecting electromagnetism. Although introducing electromagnetism 
into our calculation is not entirely straightforward because of the effects of Aharonov-Bohm phases discussed in 
ref. [T^ , we expect the approach of that paper to provide a natural generalization to the full electroweak interactions. 
In particular, the symmetry operator we rely on for the partial wave decomposition of our phase shifts continues to 
hold in the full theory. 

The paper is organized as follows. In Section II we will consider the Higgs-gauge sector model, consisting of a 
doublet Higgs coupled to an SU{2) gauge field, and discuss various aspects of classical string configurations. In 
Section III we discuss the coupling of the fermions to string configurations and compute the energy that arises from 
summing the fermion zero modes. We will present our numerical results for the corrected strin g en ergy in Section IV. 
Some preliminary results of the current investigation were reported in conference proceedings |2C| . 



II. THE HIGGS-GAUGE-SECTOR 



Our principal aim is to study string-like configurations within the standard electroweak SU{2) x U{1) model. In the 
present paper, we will consider a simplified version of this theory, in which 

1. the spacetime dimension is reduced to Z? = 2 + 1 or, equivalently, bosonic configurations are translationally 
invariant in the third space direction, 

2. the Weinberg angle, is set to zero, i.e. the U{1) factor in the gauge group decouples (and will be discarded 
in the following). The classical energy actually decreases when 6'w is turned on in the full theory. Hence a 
configuration that is stable at (?w = most likely is also stable otherwise. And 

3. quantum corrections to the bosonic string configurations are taken from the fermion sector only; this approxi- 
mation may be justified in situations where the degeneracy of the fermion modes (e.g. the number of colors Nc) 
is large. 

The action of the simplified model consists of two parts, S = Sh + Sp, where the bosonic (Higgs-gauge) sector Sh is 
treated classically and the fermion part gives rise to quantum fluctuations around bosonic field configurations. In the 
present section, we discuss properties of classical string-like configurations in the bosonic sector while the quantum 
treatment of the fermions is deferred to section IIIII 
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The Higgs-gauge sector Sh in our model reads 



(2.1) 



We work in Z) = 2 + 1 spacetime with Minkowski signature and include the gauge coupling strength, explicitly in 
all vertices. Thus, the field strength tensor is given by 

G^. = - d.W^, ~ ig [W^^, W,] , (2.2) 

where the three SU{2) vector bosons, W^, are components of a Lie-algebra valued matrix connection, 

M^^ = w^;y = w;t°. (2.3) 

As indicated, we use hermitian generators T'^ — t"" /2 (where are the Pauli matrices) with the commutation relation 
[T'^,T^] = ie^^^T". The complex scalar Higgs field is in the fundamental representation of SU{2), 



(2.4) 



while the covariant derivative contains the Higgs-gauge coupling constant g explicitly, 

D,,(b^id^,-igWf,)<j). (2.5) 

By construction, the action eq. H2.1|l is invariant under SU{2) gauge rotations V — exp{i9°' {x)T'^) , 

(j)^V(j) and igW^,^V {-d^,+igWt,)V'^ . (2.6) 

For positive values of the self-coupling A > 0, the Higgs doublet acquires a vacuum expectation value (VEV) conven- 
tionally chosen as^ 

{4>^)=v, (<^+)=o. 

As a consequence, the system undergoes complete spontaneous symmetry breaking SU{2) {1}, i.e. there is no 
group generator that leaves the vacuum invariant. As a consequence, there are three would-be Goldstone modes from 
(f> fluctuations, which become longitudinal modes of the gauge field as it acquires a mass Mw — gv/^/2 (at tree level). 
The remaining Higgs excitations are also massive with a (tree-level) mass of Mh = 2v^/X. 

In the course of this paper it will occasionally be convenient to re-write the covariant derivative, eq. (|2.5|) as 

D^<f> = [a^ - igQzZ^ - ig {W+T' + M^~T+)] , (2.7) 

where the charge operators and fields are 

' f (2.8) 

z = w3 ^ (wi ± iw^) . 

We will use the notation Wf^ and interchangeably. 



^ In the full 5(7(2) X U{1) model, this choice ensures that the unbroken U{1) group is generated by the electric charge Q = + Yw 1'^, 
where 1^ and weak isospin and hypercharge, respectively. 
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A. The Classical String Solutions 



Our main interest is in string-like configurations of the bosonic sector. We analyze the energy of (classical) localized 
configurations in two space dimensions. A subsequent extension of this analysis would be to consider configurations in 
three dimensions that are translationally invariant in the z direction, with a finite energy per unit length in a localized 
region of the x — y plane. The region in which the energy density is concentrated would thus take the form of a tube 
extending along the z-axis. Such field configurations are commonly called cosmic strings to distinguish them from the 
fundamental objects in string theory. In the present paper, we will refer to these objects simply as strings or vortices. 

Finite energy requires the configuration to be pure gauge at spatial infinity. Let (p, if) denote the radial and angular 
coordinates in the plane. At p ^ oo, the most general vacuum configuration is 

and gWjf^^ ^iVd^V'^ , (2.9) 

where the SU{2) rotation V{ip) maps the circle Si at infinity onto the gauge group G. For the electromagnetic case 
G ~ ^^(1), such maps would be characterized by an integral winding number ni(L/(l)) = TL and we would obtain 
topologically stable Nielsen-Olesen vortices |23|. However, the weak isospin group in our model is G = SU (2). We can 
still embed Nielsen-Olesen type solutions in this larger gauge group by considering maps from the spatial boundary 
S\ to various J7(l) subgroups of G = SU{2). There is, of course, considerable arbitrariness in this procedure, since 
no U{1) subgroup survives the electroweak symmetry breaking, so there is no 'natural' candidate for the target C/(l). 
Moreover, from Ili{SU{2)) — 0, no such embedded string configuration is topologically stable against deformations in 
the full SU{2) group; this will be discussed in more detail in section III Bl 

First we consider the U(l) subgroup generated by the charge Qz associated with the neutral (or Z^) gauge 
fields, as defined in eq. 12.8|l . In the full SU{2) x U{1) electroweak theory, Qz becomes the charge of .Z^, the gauge 
field that couples to the weak neutral current, and the vortex associated with Qz becomes the Z-string. The maps 
that specify the asymptotic pure-gauge configurations are 

y(^) = e"^™'^'^^ , (2.10) 

where n must be integral to obtain single-valued Higgs fields. To construct the corresponding vortex configuration, 
we first adopt the Weyl gauge 

W^o" = . 

In non-Abelian gauge theories, there might be topological obstructions to implementing this gauge [2^ |2^: such 
defects are, however, associated with non-trivial (periodic) boundary conditions in the time direction. For static 
configurations, no such obstructions exists and the Weyl gauge may always be attained; the residual gauge freedom 
consists of all time-independent gauge rotations. We shall adopt this gauge throughout the rest of the paper. 

The explicit form of the /Z-string now dwells in the neutral component of the gauge field {(p is the unit vector 
in azimuthal direction) 

2n 

<l^o = vfHip)e''''' and ^ — fcip)^ , (2.11) 

9P 

with all other fields vanishing. The radial functions fn and fa must vanish at the origin to avoid singularities, and 
approach unity at infinity to have finite energy. Substituting this ansatz into the action, eq. (|2.1|) , yields the classical 
energy 



Ec\ = 27r / dpp 







{^) + Wnf + -.fh (1 - faf + Xv^ (1 - fnY 
9 \ P J P 



(2.12) 



Minimizing this energy functional leads to the Nielsen-Olesen differential equations 

1 T7^ 

^ = fH + -Ih - -Ml - IgT + 2Xv'f„{l - fl) , 
P P 

- f^--f'G + li9'+9"yfUl-fG), (2.13) 
P ^ 

which can be integrated numerically. Near the origin p — s- 0, the radial functions behave as 

fH{p)<^P and /G(p)ocp2. (2.14) 
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Notice finally that our / Z-voitex carries quantized magnetic flux associated with the Abelian charge Qz- Defining 
the pseudoscalar magnetic field B = dxWy — dyW!^, we have 

F= fd':,B = ^. (2.15) 
J 9 

This flux would flow along the vortex tube in Z? = 3 + 1 dimensions. The classical energy density creates a pressure 
that tends to spread out the flux, while the Higgs condensate is suppressed in the region of non- vanishing _B, which 
tends to compress the flux. As a result of these two competing effects, the / Z -sirmg is stabilized at an equilibrium 
thickness. 

In the above exploration we have chosen to embed the string in the / Z-U{1) submanifold. This choice, however, 
is arbitrary. For example, we could equally well chose a family of U{1) submanifolds that is defined by the generator 
h ■ T where n is a unit vector. 

In addition to strings carried by the neutral gauge field, there are also string-like configurations that live entirely 
in the charged sector with = 0. 



B. Connecting Strings to Vacua: The Sphaleron Square 

String solutions are characterized by the homotopy group Yii{G) of their boundary condition H2.9|l at spatial infinity. 
Since Iii{SU{2)) = 0, no topologically stable string solutions exists in our model (nor in the full electroweak theory). 
In other words, the string solution can be continuously deformed into the vacuum with finite energies along the entire 
path of deformation. In the following we will consider such deformation paths, because they represent the starting 
point for a variational approach to minimize the string energy when the vacuum polarization from the fermion sector 
is included. 



1. String deformations 

We construct non-contractible two-parameter loops of configurations with W-strings as the top of the tightest 
loop Let /?2, ck denote angular coordinates which describe a three-dimensional sphere, ^3, with < /3i < tt 

and < a < 27r. The sphere 5*3 can be embedded as a unit sphere in a four-dimensional Euclidean space and is 
described by unit vectors 



/ sin pi sin P2 cos a ' 
cos Pi 
sin Pi cos P2 
\ sin Pi sin P2 sin a , 



(2.16) 



First, a map with unit winding number from ^3 to SU{2) is given by 

U'^^\Pi,P2,a) ^h^l-in-f . (2.17) 

To construct a map with winding number n, we may either raise the above C/^^-* to the power of n or, even simpler, 
scale the azimuthal angle a of C/^^^ ,26.] . 

{pi ,p2,a) = C/(i) {Pi,P2,na) . (2.18) 

If we now identify the subspace spanned by a with the circle at infinity in the spatial x — y plane (i.e. we set 
ip = a), we obtain a two-parameter family of maps characterizing possible string boundary conditions 1)2. 9|l : 

%,«.(^) = f^^"H6,6,'^). (2.19) 

As before, we construct string configurations with these boundary conditions by dressing the asymptotic form, eq. (|2.9|l . 
with appropriate radial functions, 

^^fHip)U^,,Uf){l) and W^^^fG{p)U^,,iA^)d^Ul^^M^ (2-20) 
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The radial functions go from zero at the origin (to ensure smoothness) and to unity at infinity (to ensure finite energy) . 
Substituting the exphcit expression for U^-^^^^, we get 



W+ = 



2n 



'i cos ^1 — sin cos ^2 
sin 6 sinC2e"^ 



9P 
V2n 



9P 



Vfcip) sin^Ci sin^6 , 

/g(p) sin^i sin ^2 (icos^i + sin^i cos 6) 



(2.21) 



For ^1 = ^2 = 7r/2, we have the W^^/Z-string configuration with the radial functions satisfying the Nielsen- Olesen 
differential equations. On the boundary of the square spanned by ^1 and ^2, all fields are zero except for the charged 
scalar which takes the form: 



6 : (t>+= ivfHip) , 



= TT : (/)+ = -ivfHip) > 



6 = 



IV e 



Mp) , 



6 



IV e 



fHip). 



(2.22) 



To summarize, we have a two-parameter family of finite-energy string configurations that we will refer to as the 
sphaleron square. The center of the square at (^ = ^2 — '"'/S) corresponds to the W^^/Z-string discussed earlier. 
From the /Z-stiing we may go to any point on the boundary of the square (say ^ = 0) and then deform fnip) 
smoothly to unity everywhere so as to finally reach a vacuum configuration. The energy remains finite during the 
entire deformation. In fact, the energy density on the sphaleron square is 



P 



sin2 ^1 sin2 ^2 + {f'H) + ^v' (l - fn) 



(2.23) 



which is maximal at the / Z-stfing, thereby justifying the "sphaleron" designation. Note that it is possible to 
extend this construction to the full electroweak theory; this approach will be used in a forthcoming study on the 
Z-string in 3+1 dimensions. 



2. Charged scalar condensate 



We make one simple expansion of the sphaleron square ansatz by introducing an additional radial function for the 
upper component of the scalar field. We will observe below that such an extension is necessary to avoid technical 
problems when integrating the second order Dirac equation for fermions in the string background, cf. sec. IIII Fl In 
addition, this modification allows us to probe the stability with respect to variations in the upper Higgs component. 
So our full ansatz for the string background in the bosonic sector reads: 

fH{p){~i cos ^1 - sin 6 cos 6) + fp{p) 
fHip) sin a sin ^26™'^ 



2n 

— <^/g(p) sin^ a sin^ ^ , 
9P 



V2r. 
9P 



-ipe 



-imp 



fcip) sin a sin 6 (icos^ + sin^ cos a) , 



(2.24) 



where fp{p) is non-zero at the origin and goes to zero as p ^ 00. The magnitude of the Higgs field takes a simple 
form along the Hne ^2 — tt/2, 

I'/'l' = Jh + fl ■ (2.25) 
In the following, we shall restrict our variational search to ^ = 7r/2. The additional Higgs component then contributes 



'fl [fl + 2fl - 2] 



(2.26) 



to the energy density, eq. (|2.23|) . 
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III. FERMIONS ON STRINGS 



As we have seen in the last section, electroweak strings are classicahy unstable. It is conceivable, however, that their 
quantum interactions with leptons and quarks might lead to stabilization: Since the Higgs condensate is suppressed 
in the vortex core, the fermions become effectively massless on the string, so energy can be gained by populating these 
states. If the energy gain is large enough, the string ceases to decay and it becomes a stable multi-quark object with 
a potentially rich phenomenology |^ |y| ■ 

To treat the fermions consistently, we not only have to consider the bound states but also the full Dirac sea, i.e. 
the fermion vacuum polarization, iJvac- This quantity is the most cumbersome part of the computation and to our 
knowledge, it has not been calculated properly before. The main difficulty arises because the (string) background is 
of non-perturbative nature while the renormalization conditions are formulated within perturbation theory. In this 
chapter we will discuss our approach to this issue in greater detail. 

Let us assume that Nf fermions (with perturbative mass ruf) are trapped along the string by occupying Nf of the 
bound states that are induced by the string background. These bound states are characterized by energy eigenvalues^ 
< Ci < nif. The explicit occupation of Nf bound states then contributes 

Nf 

i?^^/'=E^» (3-1) 

1=1 

to the total energy. It is, of course, clear that the minimal total energy is found if we occupy the Nf lowest bound 
states in eq. The total energy associated with fermion number Nf is then 

Egg^'' = Eel + Eo'^J^ + E^s.c ■ (3-2) 

We emphasize that Eo^J^ and -Evac are formally of the same order in and the inclusion of either thus requires the 
other. The effective energy £^g^^^ must be compared with the mass of Nf elementary (perturbative) fermions, i.e. 
stabilization occurs if 

E'-^'^KrufNf. (3.3) 

Our main goal is to perform a variational search in our string ansatz space, so as to minimize until eq. 

holds. It is immediately evident that certain situations will be favorable a priori: 

1. The perturbative fermion should be as heavy as possible; we will therefore concentrate on the top quark contri- 
bution to E^^'\ 

2. If fermions carry an internal degree of freedom (e.g. color) and Nc is the corresponding degeneracy, then 

both E'occ/'' and E'vac acquire a factor Nc, whereas Ec\ does not. Thus a large color degeneracy enhances the 
quantum piece of the total energy. The fermion number with internal degeneracy is then (Nf Nc), and the 

stability criterion eq. H3.3|l becomes E^g^'' < ruf [Nf Nc). 



A. The Theory 

We ignore inter-generation mixing and set the CKM matrix to unity. For simplicity, we shall also assume that 
fermions within an isospin doublet are degenerate. This simplifying assumption is violated in reality (in particular for 
the heavy quarks), but it gives us an additional symmetry which we need to make the calculation of i?vac manageable. 
Let us therefore consider one degenerate isospin doublet \1/ of heavy fermions (e.^. the top-bottom pair). In D = 2 + 1 
the action for the doublet is 

W^,] = / (fx [mYD^PL-^ + mj^df^PR^ - {^Pr + , (3.4) 



^ Our background profiles lead to CP invariant interactions for the fermions, cf. sec. 1111 It is therefore sufficient to consider non— negative 
energy eigenvalues and positive fermion numbers Nf. 
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where 

It 7' 



L,R 



and !t . (3.5) 



The degeneracy of the isospin doublet is reflected by the appearance of a single Yukawa coupling constant, /. The 
corresponding fermion mass, 

ruf = fv (3.6) 

is generated as usual from spontaneous symmetry breaking via the Higgs VEV. The coupling to the gauge fields 
appears through the covariant derivative, which takes the same form as eq. (j2.5(l . except that the gauge fields only 
couple to the left-handed fermions Pl"^ (as in the full electroweak theory). Finally, we employ 4-component Dirac 
spinors even in D = 2 + 1 to mimic the full D — 3 + 1 calculation as closely as possible. Including the two isospin 
degrees of freedom, our fermion spinors ^ thus have 8 complex components. 



B. The Effective Energy 

In path integral language, the vacuum polarization energy is obtained by integrating out the heavy fermion doublet. 
For this purpose, we decompose the fermion action l|3.4(l into the free part and the interaction with the bosonic 
background fields. The latter enters in the form of a background potential, U = U{W, $), 

S'f[«', Wf,, = y" [^(n^^M - rufl)^ - ^ U{W, $) , (3.7) 

which may easily be extracted from eqs. H3.4|l and (|3.5I) . If we normalize the fermion loop contribution to the boson 
action by subtracting off the result in the free non-interacting case, a formal path integral expression would be 

SM, ^] - SM = (-.) log M^^^UM^sA^,u]) 

' ^ 9L J V ; J di'f,^] expiiSpi'f^U = 0]) ^ ' 

In the language of Feynman diagrams, eq. H3.8|l equals the sum of all graphs with one fermion loop and an arbitrary 
number of {—ill) insertions, i.e. 



SM^Y.^'fM- (3.9) 

n=l 



Notice that the n = term (corresponding to a cosmological constant) is not included due to our normalization 
eq. H3.8|l . As it stands, eq. (|3.9|l is, of course, ill-defined due to severe UV divergences in the low-order Feynman 
diagrams. We can identify these divergent graphs by power counting and assign a definite value to them using a 
specific regularization method that introduces a regulator (for instance a momentum cut-off. A, or the deviation, e, 
from the physical dimension). Next, we cancel the divergences by adding countcrtcrms to the initial (tree- level) boson 
Lagrangian. If the regularization method respects the gauge invariance of the fermion loop, the counterterms will 
be gauge invariant, too.'^ Moreover, since our model is renormalizable, the counterterms are local monomials in the 
boson fields of mass dimension four or less. With these constraints, the most general counterterm Lagrangian in the 
Higgs-gauge sector reads 

ifi;*^ = citr (G^'-G,,) + C2 + C3 [{^U) - v^] + c, [(0+0) ^ y^]\ (3.IO) 

After adding this term to the initial Lagrangian eq. (|2.1|l and replacing all fields and couplings by the renormalized 

ones, the fermion loop calculation is now based on the starting bosonic Lagrangian ifl'^^+^if. The counterterm 
coefficients Ci are of order 0{h) and enter the one-loop effective action directly, 

S^r\^^^] = Sh[W,^] + ( 4'*^[V7,$] + 4d[C^] ) = SH[W,<f] + 5vac[W^,*] . (3.11) 



^ This statement is more subtle when boson loops are included. In this case, the gauge fixing requirement only leaves a global BRS 
invariance of the bare action, which does not translate directly into BRS invariant counterterms. 

If the context does not lead to confusion, we will use the same notation for renormalized and bare fields/parameters. 
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We can adjust the coefficients Ci sucli that they cancel the divergences in the low order Feynman diagrams. The 
remaining (finite) arbitrariness is fixed by imposing specific renormalization conditions on Green's functions at pre- 
scribed external momenta (discussed in more detail in section Fill D|l . Although there is no unique specification of the 
theory by renormalized parameters, we prefer to choose a scheme where these quantities take empirical values that 
may, in principle, determined from experiment; a complete list is given in the next subsection. 

The discussion above concentrated on the effective action. For static boson fields, all these statements translate 
rather trivially into expressions for the energy: We simply have to replace S —^ E = — limr^cc S/T, where T is the 
time extension of our fields. For instance, the one-loop effective action gives rise to a one-loop quantum energy 

E,s[W, $] = E,,[W, + E^,,[W, = - lim T'' S^'^''\w, $] (3.12) 

T — >oo 

where E^i refers to the classical energy of the Higgs-gauge sector, cf. eq. 12.23|l . The vacuum polarization piece in 
eq. H3.12|l has two parts, the total fermion loop correction (|3.9|l (with the replacement S E ^ —S/T) and the 
counterterm contribution 

E,t[W,^] = J d2^|-citr(TyyWy) + C2[A0]^A0-C3 [{cb^) - v^] - C4 [(0^) • (3.13) 

It should be emphasized again that only the sum of these two parts is finite as the regulator is removed, e.g. A — > 00 
or e — !■ 0. One of the main advantages of our method is that the two contributions to -Evac are manifestly combined 
such that no explicit cutoff is required even at intermediate stages of the calculation. 

Besides the classical, vacuum polarization and counterterm piece, a given configuration with fermion number Nf 
will also have to include a contribution from explicitly occupied bound states, 

E^^^f^W, $] = E,i[W, $] + E^^J\W, $] + Svac[T^, $] . (3.14) 

This is the quantity to be used for the main stability criterion eq. (|3.3|) . In the following subsections, we will discuss 
the individual ingredients into our main formula (|3.14() in more detail. 

To close this section, a brief comment is in order: We are well aware that a consistent one-loop calculation must 
include gauge and Higgs field loops, too. However, we choose to ignore these contributions in the present calcula- 
tion, since our main motivation is slightly different: We are mainly concerned with the study of fermion effects on 
bosonic backgrounds. A consistent treatment of the fermion sector requires the inclusion of both bound state and 
sea quark contributions. Although bosonic loops also appear at the same order in h, they are suppressed compared 
to the fermionic fiuctuations by a factor of the number Nc of internal fermion degrees of freedom {e.g. color). Our 
approximation is relevant when Nc can be assumed to be large, becoming exact in the limit Nc — > 00, since all boson 
lines in Feynman diagrams are suppressed in this case. 



C. Theory Parameters 

As mentioned earlier, we choose to define our model by renormalized parameters that have, in principle, a direct 
relation to electroweak phenomenology. However, our model is defined in two space dimensions so that the (dimen- 
sionful) gauge coupling constant g can no longer be expressed in terms of the Fermi constant Gp ~ 10~^m~^^^^^ via 

the standard relation Gf/\/2 = g"^ /M^ = 2/u^. To see which value of v should be chosen in D = 2 -f 1, consider the 
classical energy (or energy per unit length) in eq. (|2.23|l . In any number of dimension this quantity takes the form 
times a dimensionless function involving only ratios of masses. Thus the mass dimension comes from the prefactor 
z)^ only, which must therefore be adjusted such that the known D = 3 -I- 1 value of v'^ is multiplied by an appropriate 
length scale. Since we relate all dimensionful quantities to the mass mj of the heavy fermion, one might expect that 
the appropriate length scale is the Compton wave length of the heavy fermion, 27r/m/. This is almost correct; the 
actual calculations of quantum energies for the QED fiux tube (eqs. (42) and (43) in ref. show that the proper 
scaling factor is in fact half the Compton wave-length. Thus, the physically sensible choice for is 

2 TT 2 

^D=2+i — ^D=3+i = TTIkT^ ■ (3.15) 

TO/ 2V2GFmf 

For the remaining parameters, i.e. the masses that determine the coupling constants /, g and A, we can directly take 
the values suggested by electroweak phenomenology: 

m/ = /t)i)=3+i = 170GeV, = 80 GcV and 2wd=3+i\/A = 115 GcV . (3.16) 

v2 

We have identified the heavy fermion doublet with the top/bottom pair and set mj to the top quark mass. In our 
numerical studies, we have also experimented with a very heavy fermion to investigate possible decoupling scenarios. 
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D. The Counterterms 



The counterterm coefficients Ci in (|3.13() must be chosen to cancel the divergences in low-order Feynman diagrams. 
In D = 2 + 1, only the first two diagrams (with n = 1 and n = 2 insertions of the background potential —iU) are 
divergent and it is readily seen that the terms in (j3.13(l are sufficient to render the sum of the n = 1,2 diagrams — 
and thus the entire vacuum polarization part i?vac — finite. The remaining arbitrariness in the definition of the c; is 
fixed by conventional renormalization conditions (see, for instance, ref. 27]): 

a. We start by choosing the so-called no-tadpole condition. This ensures that neither the VEV (0) — v {0 , 1)^ 
nor the perturbative fermion mass ruf receive radiative corrections. This condition determines C3. 

b. We fix the pole of the Higgs propagator to be at its tree level mass, Mh, with residue one. These conditions 
yield C2 and C4. 

c. There are various choices to fix the remaining coefficient ci. We choose to set the residue of the pole of the 
gauge field propagator to one in unitary gauge. The position of that pole, i.e. the mass of the gauge field, Mw, 
is then a prediction. 

The resulting counterterm coefficients Ci are listed in the Appendix. As explained under item c, the mass of the 
gauge fields is constrained by the other model parameters when fermion loops are included. With our choice of 
renormalization conditions, Mw is the solution to the implicit equation 

= S(l + -^HCw)) + T/^(e^/i(ew) + 2/1 (^w)), (3.17) 
2m j: nmf loirmf 

where = Afw /"tl / ! the Feynman parameter integrals Ii , Ii and I2 are also given in the Appendix. If we take for g 
and mf = fv the values extracted from eqs. (|3.15|l and H3.16|l , we can solve eq. (|3.17() numerically to obtain the fermion 
loop correction to the gauge boson mass. Due to the condition c. above, this mass corresponds to one-Vt^-boson states 
in unitary gauge. 

Of course, other renormalization conditions are possible. They correspond to finite changes in the parameters 
a which are exactly calculable and affect the counterterm contribution eq. (|3.13|l only. The renormalization group 
ensures that such re-parameterizations do not change the physical content of the model. 



E. The Vacuum Polarization Energy 

From eq. H3.9|l the one-fermion-loop correction to the bosonic background energy is an infinite series of Feynman 
diagrams, which is impossible to sum in practice. Perturbative methods based on the lowest order diagrams are 
ineffective, since the coupling of the heavy fermion to the boson background is not naturally small. Instead, we 
follow the spectral approach introduced in refs. '2^, '29'| for an exact calculation of the renormalized fermion vacuum 
polarization. For a detailed review of this method and a list of further references, see ref. |l_7| . 

The spectral approach is based on the observation that the vacuum polarization can alternatively be characterized 
as a sum over the change in the zero-point energies of the fermion modes due to the background fields. The fermion 
modes comprise both the bound states of energy ej (with degeneracy Dj), and the scattering modes characterized 
by a momentum density of states. If Ap(fc) denotes the change in the density of the scattering states caused by the 
non-trivial background, we have the expression 

1 1 f°° I 

where the degeneracy factor Dj m D = 2 -\- 1 dimensions is unity if the sum over j runs from —00 to 00. The 
cancellation of the divergences in the momentum integral by conventional counterterms i?ct is still formal at this 
stage. 

What makes the spectral method effective is the fundamental relation between the change in the continuum density 
of states and the scattering phase shifts 5(,{k) of the Dirac wave-functions [s^ 

Mfc) = ^^E^^^^W- (3.19) 
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The shorthand notation £ wiU be used in the foUowing to denote the coUection of aU quantum numbers (including the 
sign of the energy and generahzed angular momentum) that characterize an individual scattering channel. 

Since the Feynman series is an expansion in powers of the interaction potential eq. (|3.7(l . the same expansion in 
the spectral method suggests a one-to-one correspondence between Feynman diagrams and the Born expansion of 
the phase shifts. This reasoning is essentially correct, though there is a slight subtlety associated with the string 
background. Ignoring this issue for a moment, we may subtract the lowest orders of the Born series from the phase 
shifts, and add back in exactly the same quantity in the form of Feynman diagrams: 



1 



2 ^ 



D,{\e,\ - 



rrif - - 



N 



(3.20) 



Here, [6{k)] ^ is a shorthand notation for the sum over all channels £ of the full phase shift (includin g the degeneracy) 
minus the first N terms of the respective Born series. Ref. ,31j contains a careful derivation of eq. (|3.20l) that starts 
directly from the definition of the energy-momentum tensor in field theory. 

All three contributions in eq. H3.20|l are separately finite as the regulator is removed: The bound state is a finite 
sum, the momentum integral converges due to the Born subtraction, and the last piece is a conventional perturbative 
calculation of the N lowest order Feynman diagrams (with counterterms to implement the renormalization conditions 
discussed above). For the present model in D = 2 -I- 1, only the first and second order Feynman diagrams require 
renormalization, i.e. N — 2 is the minimal number of Born subtractions necessary for a finite calculation. The final 
expression for the vacuum polarization energy is thus 



E\rp,c 



1 , ^ dk , [ 



k^ + mj-mf)— [S{k)]2 + E'^'^'' + E, 



p(3,4) 



(3.21) 



The first and second order diagrams, combined with the counterterms of the same order give the finite contribution 
^(ii2)_ fpj-^g third and fourth order diagrams are not divergent in _D = 2 + 1 (and therefore left as implicit part of 
the momentum integral), but we still need finite counterterms £'ct'^'' for them to enforce the correct renormalization 
conditions. The explicit expressions for E^^''^^ and -E^t can be found in the Appendix. 

Let us finally discuss the subtlety associated with the Born approximation in the string background. It is due 
to the fact that Feynman diagram calculations require all field derivatives to vanish sufficiently fast as p ^ oo. 
This condition is not met by the vortex, which has winding and thus nonzero azimuthal derivatives at p cxd. In 
D = 2 -|- 1, only the counterterm coefficient C3 is UV divergent while the remaining (finite) Ci serve to enforce our 
renormalization scheme. The counterterm multiplied by C3 involves only (f>), so instead of doing the full Born 
subtraction and the corresponding diagrams, we can perform this piece of the calculation in a fictitious background 
with the same value {4>^4>) as our initial vortex background. On the sphaleron square, we can, for example, choose 
the corresponding boundary point with ^1 = and ^2 = 7r/2, cf. Sec. Ill Bl where all gauge fields vanish and cj) has no 
azimuthal dependence at p — )■ 00. For this fictitious background, the subtraction/diagram part of the calculation is 
straightforward. 



F. The Dirac Equation 

The final ingredients in our computational method are the equations of motion for the fermion fields and the associated 
techniques used to extract the phase shifts in the boson background. 



1. First- order equations 
We start from the Dirac equation in 2 -I- 1 dimensions, which we write in the form^ 

h-D^oj.Ne (P, </5) = l^'^uj,Ni (p, f) , 



(3.22) 



^ With four-component Dirac spinors, the same equation also holds in D = 3 -I- 1, provided that all fields are translationally invariant 
along the string (which we put along the 2;-axis) and the gauge fields have vanishing z-components. This means that the bound states 
and phase shifts co mpu ted in this section are also appropriate for _D = 3 -I- 1 and the trivial space dimension can be managed by the 
interface formalism Il9l . Complications arise in D = 3 -f 1 only because the renormalization procedure is more elaborate. 
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The single particle Dirac operator can be read off from eq. H3.4|l , 



ga 



Z,Qz + Pl{W+T- + W-T+) 



(3.23) 



where a* — 7^7* and i ~ 1,2. Our vortex configurations are special cases of a class of background fields which show 
an angular dependence similar to the Nielson-Oleson vortex of winding number n, 



cxe"™"^^, 



(3.24) 



&o cx e""^ , 0+ oc 1 , oc ^ , W+ = {w~^ ' 

This class of configurations exhibits a symmetry: Since the heavy fermion doublet is degenerate, the operator 

K^L^ + S^- nT3j5 (3.25) 
2- may therefore be used to label the fermion modes 



commutes with the Hamiltonian hj). Its eigenvalues Ni 
and channels in a vortex background of the specified form. 

As explained earlier, we want to keep our computation as close as possible to the D = 3 + 1 case, and therefore 
employ Dirac 4-spinors even in D — 2 + 1. This redundancy implies the existence of a matrix with the property 
{a^, hu} = 0. For every eigenspinor 5* of the Dirac operator with eigenvalue w, there exists an corresponding state 
a^vff -vvith eigenvalue —to, and the spectrum of ho is manifestly symmetric about zero. 

For the actual calculation we choose the chiral representation of the Dirac matrices (with 7^ — diag(— 1, —1, 1, 1) 
and 5*3 — diag(l, —1, 1, — 1)/2) and parameterize the eigenspinors by eight functions yi, . . . , ys that only depend on 
the radial coordinate p. 



uj,Ni 



iy^ e*(^+i)^ 1 1) + iy^ e'^t+n+iW | _ i) 
V iy2 e»(^+"+i)v' \\) +iyi e*(^+i)^ | - i) / 



(3.26) 



The kets in this expression are eigenvectors of T'^ = /2. We may re-cast the D = 2 + \ Dirac Hamiltonian in 2 x 2 
block form as 



hii hi2 
^21 ^22 



(3.27) 



In the chiral representation of the Dirac matrices (and in temporal gauge) these 2x2 blocks have a very simple 
dependence on the background fields. 



(3.28) 



/ill 


= icTidi - 






= -icFid 




hi2 






h2i 


= /^^ 





2. Second-order Equations 



To extract the scattering data, the second order form of the equations of motion is more appropriate. We find these 
by substituting the Dirac operator eq. H3.27|l into the first order form eq. H3.22|l which eliminates half of the field 
components (while doubling the order of the differential equation). We choose to eliminate the left-handed fields, 
which are the upper components in the chiral representation. The resulting equation of motion for the right-handed 
(lower) Weyl spinor '^r reads 

/i2i [(/ill - uj)h^i {h22 -uj)- /112] ^-fl = . (3.29) 

These are two complex second-order equations for the two complex components of ^! n. Upon inserting the parame- 
terization eq. H3.26|l . we obtain four second-order equations for the radial functions appearing in ^> j^, 



|f^(p) + ^(p) A + ^(p)| y^{p) ^ 



(3.30) 
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rnj. The kinetic operator ^ contains a radial piece from the 2D Laplacian, as well as a generalized 
centrifugal barrier ^, 

^(p) = i(^ + l^ + kA ^? = diB.g[ie + n)\{e + n + i)\e\ii + if] . (3.31) 

\ ap^ p dp J 

The matrices ^{p) and ^{p) contain the bosonic background potential. The explicit form of the matrix elements 
is rather lengthy, particularly for the non-trivial background ansdtze in our numerical treatment; we refrain from 
presenting them in any detail. At large radii p — > oo, the background potential vanishes and ^{p) ^ as well as 
^{p) 0, and the four differential equations decouple. 

The bound state energies ej can be computed directly from eq. 13.22|l by shooting for normalizable solutions that 
are regular at the origin and decay exponentially at large radii p — > oo. By contrast, the phase shifts, which are 
extracted from eq. I|3.30|l . require a few more manipulations to put the calculation in a manageable form. 

We have a four-channel scattering problem that can be treated in the usual fashion: The full eigenspace to a given 
energy uj is decomposed into channels labeled by the integral angular momentum I. In each such channel, we have 
four independent solutions y{p) to eq. (|3.3U|) . which we may combine in the rows of a matrix '3^. Thus, i^ijlpjk) is 
the radial function yj{p) in the i**^ linearly independent solution at a fixed momentum k and channel number £ (we 
suppress £ in the following). 

Let us analyze the free case ^ = ^ — first. The solutions with outgoing spherical wave boundary conditions at 
p — > oo, are W {p, k) = J^{kp) with 

^(x) = diag [/jW „| {x) , /i[;|„+i| (x) , h\lj (x) , h\ll^^ (x)] (3.32) 

given explicitly in terms of Hankel fimctions. Similarly, for ./#, ^ ^ 0, the scattering states y(p) are solutions 
of the fully interacting equations of motion H3.30|l . which tend to specific free outgoing waves of the same energy 
at p — > oo. Combining them in a (4 x 4) matrix as before, we may split off the asymptotic form by the ansatz 
W (p, fc) — ^(fc, p) • J^{kp). Inserting into cq. H3.30|l . we find the differential equation for the Jost-like matrix =^(fc, p), 



p 



^" + 2^'Jf' ^[JT, 0\ -f Ji{^' -f -t- = , (3.33) 



where J£[kp^ = liiJ^(kp) and primes denote derivatives with respect to the radial coordinate. If we impose the 
boundary conditions 

jr(fc,p)^l, .^'{k,p)^0 

at p — > oo, the corresponding solution matrix W (p, k) = ^{k, p) ■ J^{kp) clearly describes outgoing spherical waves. 
What is unusual about the string background is that the matrices ^ and yy may actually be complex. As a 
consequence, the conjugated matrix ^* does not describe incoming spherical waves; in fact, it is not even a solution 
of the equations of motion. Instead, the scattering states corresponding to incoming waves, which we denote by '3^ 
in the following, must be computed separately by an appropriate change in the Hankel functions — *■ Jff. The 
standard scattering analysis now proceeds as usual with iV* replaced by iV.^ 
The scattering wavefunction can now be written as 

3i,{p)=-W{p) + ^{p)-y{k). (3.34) 

Requiring that iV^dp) be regular at the origin yields a relation for the scattering matrix, 

J5^(fc) = lim^-i(p) -^(p) = limJ'/f-\kp)-^~\k,p)-:W{k,p)-Jf{kp). (3.35) 

For the spectral method (|3.21|) . we need the sum of the eigenphase shifts obtained from S^, 

6{k) = -^TTlny{k) = ^ \imTTln{.^-dk,p) ■W{k,p)) . (3.36) 



Alternatively, we could rewrite the complex (4 X 4) problem as a real (8 X 8) problem. The phase shifts — or more precisely the sum of 
eigenphaseshifts in each channel - extracted in this way are exactly twice as large as the ones from the initial (4 X 4) problem. This 
statement is obvious in the case that the matrices ^ and ,jV are real, and we have checked it numerically in other cases. 
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It should be noted that the fermion flux in the chiral basis is not orthonormal and the matrix ,5^ is hence not unitary. 
However, the physical scattering matrix can be obtained from ^5^ by a change of basis and proper normalization. Both 
manipulations do not alter the sum of eigenphase shifts, which is the relevant quantity here. 
In our numerical calculations, we do not use eq. 13.36|l directly; instead we set 

<5(fc, p) = ^Tr In (jT"! (fc, p) ■^{k.p)) , (3.37) 

such that (5(fc) — b(k,G). To compute this function, we integrate the differential equation 
9(5(fc,p) 1 



dp 2i 



■Tr 



lim b{k,p) = (3.38) 

p — ^oo 



along with the system for ^ and ^ from p = cx) to 0. This procedure results in a smooth function of both p and fc, 
which avoids ambiguous jumps of (multiples of) tt and allows to take the limit p ^ easily. 

Notice finally that channels which involve Hankel functions of order £ = are afflicted by a further numerical 
subtlety. Eq. (|3.36(l essentially requires to separate the regular and irregular pieces as p ^ 0. For ^ = this amounts 
to distinguishing a logarithmic behavior from a constant, which is difficult to handle numerically. In this case, we 
switch to a radial function analogous to S{k,p), which is defined in terms of the radial derivatives of .^{k,p) and 
J^{k,p). We have checked the numerical stabihty of this treatment thoroughly. 



G. Derivative expansions 



To check our numerical results, it is helpful to have analytic approximations to the vacuum polarization energy. The 
lowest orders of a perturbative expansion are already contained in our main formula eq. H3.21II ; to check the phase shift 
piece of the calculation, we would therefore have to compute (many) higher order diagrams, which is very cumbersome. 
A better strategy is to employ a derivative expansion which is expected to work well for reasonably broad background 
profiles. 

The leading order of the gradient expansion is the effective potential, T4ff, which contains no derivatives of the 
boson background fields. Together with gauge invariance, this restriction rules out any dependence on the gauge 
fields, so V^fi is a functional of the Higgs field alone. In = 2 + 1 dimensions, V^ff is formally given by 



/ 



d^x Vm 



-iTrln {i^~ /$5) 



(3.39) 



for a constant background $5 — ^Pl + ^^Pr, cf. eq. I|3.5|l . The functional trace is divergent and needs to be 
renormalized by adding the same counterterms as were determined in section Fill Dl To get definite expressions, we 
employ dimensional regularization and write 



^3^ dVoff[<l>] 



df 



-(-z)/Tr \d^ + f<S>,^l 



tr 



fc2 - /$5$t 



(3.40) 



The analytic continuation of the integral on the rhs to I? = 3 is finite, so that one would expect to need no renor- 
malization at all. However, the integral is manifestly convergent only in D ~ 1 and the continuation to Z) = 3 is 
unambiguous only if the divergence in the even dimension _D = 2 is canceled by counterterms. Only the sum of the 
integral and counterterm is analytic in the range D = 1 to D = 3 and we must renormalize for all D > 2. Furthermore, 
we need to impose renormalization conditions that match our previous scheme if we want to compare the results of 
the derivative expansion with the spectral method. Upon integrating eq. (|3.4()(l with respect to the Yukawa coupling 
and fixing the integration constant via the normalization T4fi[w] = 0, we obtain the result 



TT 



C4 



(3.41) 



The explicit expression for the counterterm coefficient 04 is rather lengthy; it can be read off from eq. (jA2|) in the 
appendix. Notice that the effective potential is minimized for configurations with (j)^(f> — u^; this is, of course, a direct 
consequence of our no-tadpole renormalization condition. 

For the next-to-leading order in the derivative expansion, there are well-known techniques [s^ Issf that apply 
straightforwardly in the case of a pure Higgs background. On the sphaleron square, these are the configurations with 
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^1 — 0. We will see in the next section that such backgrounds play a prominent role. The result for the effective 
Lagrangian is 



(3.42) 



where the superscript indicates that this term contains two derivatives of the Higgs background. Again, the expression 
eq. (I3.42|) is ultra-violet finite in D = 2 + 1, but we need to add the two derivative counterterms, eq. (|3.10|l . to match 
our previous renormalization conditions. 

In a gauge invariant expansion, the leading order effect of the gauge bosons is to turn the ordinary derivatives in 
eq. H3.42|l into covariant derivatives of the form eq. (|2.5|l .'^ Inserting the explicit form eq. 12.21|l of the static vortex 
background, our final expression for the derivative expansion of the vacuum polarization energy is 
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{fH + fpy'^'ifHf'H + fHf'Hf 



(3.43) 



where primes denote derivatives with respect to the dimensionless variable x = mfp. 



IV. NUMERICAL STUDIES 



Our goal is to find a spatially varying configuration of the boson fields that is energetically favored over the trivial 
vacuum configuration with the same fermion number. To this end we scan the energy surface with respect to variational 
parameters in specific ansdtze for the background profiles. The comparison requires, of course, that both the trivial 
and non-trivial configuration have identical fermion numbers. As discussed at the beginning of section ITTll the vortex 
background acquires its fermion number by occupying the Nf lowest fermion bound states; its energy must then be 
compared to the trivial configuration with energy Nfirif, cf. eq. H3.3|l . These considerations become slightly more 
complicated if internal degrees of freedom {e.g. color) are included, as is discussed in more detail below. 

Our variational ansatz is discussed in detail in the next subsection. Altogether, it involves five parameters; even 
though this number is already numerically expensive to be scanned, it is definitely not exhaustive, so stable configura- 
tions can still exist that we have not found. However, if we find an energetically favored object, expanding the ansatz 
will only improve on this result. We will begin by exploring the parameter dependence of the various contributions 
to the total energy and then turn to the search for bound objects with large fermion numbers. 



A. Explicit profiles 

Unless stated otherwise, we shall employ the following ansdtze for the radial profiles in the string background: 

fH = l-e-P/'"", /g = 1 - and /p = ape"''/'"^ . (4.1) 

The constants wh, wq and ap, together with ^i, eq. I|2.24|l . are variational parameters that characterize the back- 
ground fields, which we will adjust to minimize the energy functional E^g^\ eq. H3.2|l . The profile functions fn and 
fa defined above are constructed such that the classical energy is always finite, regardless of the values chosen for the 
variational parameters. Near the origin, these radial functions behave as fn oc p and fc oc while for large p both 
profile functions tend to unity. We certainly could introduce additional degrees of freedom subject to the boundary 
conditions. Introducing additional parameters for the amplitudes of the exponential functions in fn and fc violates 
this requirement of finite classical energy, however, and hence we refrain from doing so. 



There are various approaches to counting derivatives in a covariant expansion; some of these schemes expand in the gradients of the 
magnetic field, so the leading term will contain B (and hence the derivatives of the gauge field) at all orders. For our purposes, the 
simple expression i3.43i is sufficient. 
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B. Classical Energy 



We first discuss the results for the classical energy Ec\. We consider the M^-string on the sphaleron square, eq. H2.24|l . 
with the classical energy 



The classical energy density, £c\, is given by eq. (|2.23l) plus the contribution of the additional scalar field fp as given 
in eq. (|2.26l) . We always assume unit winding n — 1. 

We first discuss the case for physical (standard model) parameters augmented by the choice of the VEV in = 2 + 1 
dimensions, eq. H3.15|l . In what follows all data are measured in units of m/ or its inverse. For the classical energy, 
we observe a strong dependence on the angle ^i, which parameterizes the interpolation between the VK-string and a 
purely scalar configuration with the same Higgs profile 

Even though the mixing angle is one of our variational parameters, it is fruitful to discuss the energy functional 
for various fixed values of and attempt to minimize Ed in the space of the remaining variational parameters. Note, 
however, that minimization of Ed is not the main goal of our investigations. For small the gauge field contributions 
are suppressed and the minimal value for Ed is obtained for small Higgs field extensions. In the extreme case of no 
gauge fields (^i ~ 0) we expect Ed « 6 as wh and ap = from the scaling properties of the derivative terms 
in / pdpSd, cf. eq. (|2.23|) . For = O.OStt, wh — 1-1, wp = 2 and ap = 0.1 we find Ed ~ 9. The classical energy 
essentially increases quadratically with wh- On the other hand, the dependence on wp appears to be weak, even 
for sizable ap > 0.5. For moderate « 7r/4, the classical energy is again minimized by narrow Higgs fields while 
the gauge fields tend towards intermediate sizes {wq « 3). However, the value for the (minimal) classical energy at 
fixed ^1 ~ 7r/4 is increased by about a factor 2 as compared to ~ 0. In the regime of its maximal value ~ 7r/2, 
the (minimal) classical energy is even further increased to about 30. Again, in this case the minimal value occurs at 
non-zero extension of the Higgs field wh ~ 3. When we study the effect of the additional contribution from fp, we 
also see that the minimal classical energy occurs at small and it does not vary much in the range ap e [0.1,0.5]. 
Beyond that regime the classical energy starts rising linearly with ap. 

We find it worth mentioning that the configuration we have investigated in detail with the lowest classical energy 
does not have exactly = 0, where the gauge field part in the background configuration vanishes identically, but 
rather = O.OStt. As we have seen before, even for this small value of the gauge field energy is of a noticeable 
size. We note, however, that = is problematic for the numerical investigation as will be discussed further below, 
cf. footnote 9. 

To some extent the situation is different for large fermion masses. Even though the classical energy does not 
explicitly depend on m/, this quantity enters via the VEV in Z) = 2 + 1 dimensions, eq. H3.15|l . and also affects the 
numerical result for Ed as it sets the overall energy scale*. We have numerically studied the case with ruf = 1.5GeV. 
For small the classical energy is also small, i.e. for < 0.1 we find Ed ^ 0.4. The numerical values for the classical 
energy seem quite insensitive to the shape of the Higgs field as the coupling via the VEV is drastically decreased. Of 
course, for small the dependence on wg is again mitigated. This changes for moderate values of ^i. For « 0.5 we 
have varied wq G [1-5, 5.5]. In that range Ed strongly varies between approximately 2 and 25 with the smallest result 
obtained for the largest width wg- On the other hand, the dependences on the Higgs field variational parameters 
is soft. We observe a similar situation for in the vicinity of 7r/2, when the relative contribution from the gauge 
fields is maximal. However, the range along which Ed changes is even bigger. For small wg 1-5), it may be as big 
as 50 while it reaches values close to one for very large widths of the gauge fields wg ~ 12. We do not thoroughly 
study such configurations because in these cases the computation of the vacuum polarization energy is numerically 
very costly. Minimization of the classical energy requires a suppression of the gauge boson contribution (small ^i) 
and/or a sizable extension of these fields. In the following sections we will study the behavior of the fermion vacuum 
energy under these conditions. In particular, we will study the question of whether or not the energy gain due to the 
emergence of a fermion zero mode at — 7r/2 can compensate for the increase of the classical energy. The energy 
gain per occupied zero mode is unity and thus we expect that we need a large number of fermions (twenty or so) to 
break even, unless something unexpected happens to the vacuum polarization energy. 

As mentioned above, the classical energy does not vanish for the limiting case of our variational parameters in 
which ^1, Wh and ap are all zero. The corresponding classical energy Ed ^ 6 is the minimal value that can be 



If one re-expresses the classical energy by replacing couplings by physical masses and the Fermi coupling and takes into account that 
the energy is measured in units oi mf it is obvious that the Higgs field contribution to the energy is suppressed relative to gauge field 
contribution when nrif gets large. 




(4.2) 
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accommodated by our parameterization of the fields. For this configuration, we expect the fermion spectrum to be 
that of the trivial vacuum, hence there appears to be a minimal energy threshold that must be overcome by finite 
size background fields. 



C. Bound States 



Central to our analysis are the bound states of the W- and Z-strings. For Ci = f 0+ = 0, an exact zero 
mode exists. Contrary to early expectations, it does not stabilize the 2^-string configuration against 0+ condensation 
once effects of the fermion determinant are taken into account, as was shown in the perturbative analysis of By 
taking the whole bound state spectrum into account we find something unexpected. In the background of rather 
wide strings (several Compton wave lengths of the fermion) we find a very large number of bound states. Because of 
this large number of bound states (both spread over a lot of angular momentum channels and in individual channels) 
the properties of individual states like the zero mode are a lot less important than previously thought. This result 
provides an important reason for why the purely scalar configurations are favored over the configurations containing 
gauge fields. For the purely scalar configurations the lowest-lying states are less strongly bound, however, and 
contrary to expectations from potential theory in quantum mechanics, many more bound states exist than for those 
configurations that generate strongly bound states. Thus the decrease of binding of the lowest -lying states makes 
only a small destabilizing contribution. But the price in classical energy to pay for gauge field configurations is rather 
high, as we have already discussed. Hence it can be advantageous to switch them off. 

Before discussing the numerical results quantitatively we would like to introduce criteria for their accuracy. Nu- 
merically we find the bound states via a shooting method, in which we determine the number of bound states in each 
channel by Levinson's theorem. If our numerical precision is not high enough, the shooting method will miss bound 
states that are very close to threshold^. Since each state contributes to the total energy only through its binding 
energy, these states are relatively unimportant. Therefore it is more efficient to construct an accuracy criterion from 
that observation and reject those computations that do not satisfy that criterion. Let pi be the number of missed 
bound states and LUg be the largest bound state energy that we observed in the channel with angular momentum £. 
From eq. (|3.21|l an upper limit for the numerical error in the vacuum energy is then 

1 °° 

t=~oo 

and we reject any computation for a given set of variational parameters that exceeded A_Evac by a prescribed value, 
which we generally take as 5% of the vacuum polarization energy, defined in eq. (|3.21|l with Dj = 1. Most of the 
rejected calculations are characterized by very small ^i. 



. a) 


1 1 1 1 1 1 


Ill 




0<E< 1.0 






0.9<E<1.0 






0<E<0.6 , 
, ■' 1 


■ ■ = 8 

" 8 

, 1 
1 ° 

a 8 




1 ■' ' * 

.1 ■;. .'■ ■ 




• 1 


M' 


I! : ii : It : i: : E : 



b) 


1 1 1 1 1 1 


Ill 




0<E< 1.0 






0.9<E<1.0 


! ■ = ° - 




0<E<0.6 1 
1 1 




• 


1. ■' ' ' 

1 J s 

,11 ■' 

'1,' ... ..!. .... .... ... ... ■■! 


1, . , ■ , ■ , " 



Another problem is that a shooting algorithm might miss states in pairs if they are degenerate in energy. Such a degeneracy actually 
occurs for = because the CP symmetry is enhanced to an individual C and P symmetry. We counter this problem by avoiding the 
special case §1 = 0. 
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FIG. 1: The number of bound states as function of the width wh of the Higgs field for two different values of the mixing angle: 
^1 — 0.1 X 7r/2 (a) and = n/2 (b). We distinguish between strongly (0 < e < 0.6) and loosely (0.9 < e < 1.0) bound states. 
The entry < e < 1.0 obviously gives the total number of bound states. The towers in each entry show results that originate 
from various values of those variational parameters that are neither held fixed nor are shown along the abscissa. 

In figure ^ we display the number of bound states as a function of the width parameter wjj ■ There are n£ bound 
states in each orbital momentum channel and we display ne- First of all we see that the number of bound states 
strongly increases roughly quadratically with the width of the Higgs field. The number of bound states increases only 
gradually with increased coupling to the gauge fields as measured by , because the main reason for the appearance 
of these bound states is that a Higgs field with |$| < K^)] produces an attractive potential for the fcrmions, since 
/|$| is the effective fermion mass. The number of strongly bound states also increases only gradually with i^i, which is 
somewhat surprising since we have a zero mode for = tt/2 and ap = 0.^" However, it is not generally accompanied 
by additional strongly bound states. 

This bound state structure suggests that the regime with small coupling to the gauge bosons (<^i ^ 0) might be of 
particular interest in our analysis. We consider the bound state contribution to the vacuum energy 



2^ 

b.s. 



(my 



(4.4) 



and in figure[3we display Eb as function of wh and . We see that Eb increases with wh approximately quadratically, 
while the dependence of Eb on the other variational parameters is weaker. As an example we also show Eb as a 
function of in figure |2d) . Even though Eb is large for some of the configurations in the vicinity of ^ = ^ the 





U{nJ2) 



FIG. 2: The bound state contribution, Eb to the vacuum polarization energy as a function of the variational parameters a) 
Wh and b) that are defined in eqs. (12.241 and (14.111 . The towers in each entry show results that originate from various values 
of those variational parameters that are not shown along the abscissa. It should be emphasized here - especially with regards 
to panel b) - that our coverage of parameter space has not been uniform. For a large range of widths, e.g., computations have 
only been performed for = 0.1, 0.5, 0.9. Hence one may deduce from panel b) that the range of energies is very similar for 
^1 = 0.1 and ^1 — 0.9, but not that the range of energy values at = 0.2 is much smaller than for = 0.1 (in units of 7r/2). 

maximal values are found for small ^i. Since the classical energy grows with ^i, it is unlikely to find stable objects in 
the region of large ^i, i.e. for W dominated string configurations; rather pure Higgs configurations seem to be favored. 



D. Vacuum Polarization Energy 

We now turn to the discussion of numerical results for the vacuum polarization energy, E^nc, which we compute 
according to eq. ()3.21|l . We have accumulated data for about 1000 sets of variational parameters. Our results are 



We did not consider a configuration with = tt/2, wjj = 14.8 and small ap. Thus the wp = 14.8 entry in figure^}) misses a strongly 
bound state. 
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displayed in figure 13 for the physical value of the top quark mass mf = 170GeV. 

We note for some configurations the renormalized vacuum energy is negative and can provide some energy gain 
to the system. The maximal gain is achieved for moderate extensions of the Higgs field wh, wp '--^ 3 ... 5 while i?vac 
quickly increases when these extensions grow even further. As figureOt) indicates, the most favored value of the gauge 
boson field extension is about wq ^ 3. Both of these regions could also lead to a small classical energy. On the other 
hand we do not observe a strong dependence of i?vac on the angle cf. figure|31i)- That is, the vacuum polarization 
part of the energy is not capable of overcoming the strong drift towards the Higgs dominated configuration (^i = 0) 
caused by Ed- In addition, we observe that minimization of the vacuum energy induces a charged Higgs field, i.e. 
i^vac is minimized by ap ^ as shown in figure l^t). 

In all we find that the energy gain from i?vac is only of order —mf. The lowest value that we observed is i?vac — 
— 0.33rn/ for wh = 5.1, wp = 5.0, wa = 5.5, ap = 1.9 and — 0.1 x {tt/2), which is essentially a Higgs background 
configuration. So this energy gain alone is not sufficient to stabilize the string, but we can still populate fermion levels 
along the string to construct a stable object. 

E. Gradient Expansion 

We have already introduced the gradient (or derivative) expansion approximation as a method of checking our 
numerical results. In figure 0] we display the comparison of the full result, -Evac, which we compute according to 
eq. H3.20|l . and the gradient expansion approximation, as discussed in eq. (|3.43|l . We display these data at the same 
scale as the various pieces that contribute to E^^c to show that the sizable cancellations that occur among these 
piece are necessary to accomplish the agreement with the gradient expansion approximation. In particular the case 
ap > 0.8 shows that changes in the bound state contributions E'bs are compensated by a corresponding change in the 
(renormalized) Feynman diagram contribution, -EpD- Naively we expect compensations between the bound state and 
phase shift contributions as a reflection of Levinson's theorem. However, the Born subtractions relocate substantial 
pieces of the phase shift contribution into the Feynman diagrams. We only show the results as functions of the width 
of the Higgs field, wh, because the dependence of the subtracted phase shift contribution Es on variational parameters 
other than wh is negligible. Some small discrepancy between the -Bvac and i?DorExp occurs at small because of the 
numerical problems we encounter in finding very weakly bound states in that case. However, this is not too worrisome 
because we can easily extrapolate from the region > 0.01. 

F. Stabilization with Populating Levels 

We have seen that many fermion bound states emerge for the background string configuration. We may populate 
these levels to describe an object with fermion number NfNc and effective energy 

E^^^f^ ^NcJ2^^ + NcE^^c + E,i , (4.5) 
1=1 

where the sum over the bound states is such the Nf lowest energy modes are included. The corresponding binding 
energy is 

giNf) ^ ^(Nf) _ ^^^^^^ = NcY^ (e,; - m/) + iVc-Bvac + ^^ci • (4.6) 

1=1 

Obviously B^^f^ < indicates an energetically stable object with fermion number NfNc- If there are fewer than 
Nf bound levels for the configuration under consideration, we have to "occupy" scattering states at threshold, which 
do not contribute to the binding energy. In figure [31 we display the extremal value of the binding energy for a given 
number (Nf) of populated levels. That is, for a prescribed value of Nf we extract the lowest iJ^^/^ in the space of 
variational parameters. Bound objects exist for sufficiently large Nc- In this case Nc > 8 yields a bound object with 
Nf > 20. 

For small Nc, the effective energy is dominated by Ed, which in turn is minimized by small wh- In that case 
there are only few bound states to be populated and B^^i^ is saturated once all of them are occupied. For larger Nc, 
the fermion contribution to the effective energy becomes more important and the minimal value is obtained at larger 
widths, where there are more bound states and the saturation of sets in at larger values of Nf. This behavior 

is also refiected in figure |S| 
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FIG. 3: The vacuum polarization energy as a func- 
tion of the variational parameters a) wh, b) wp, c) 
WG, d) ^1 and e) ap that are defined in eqs. 12.241 
and H4.1|l . The towers in each entry show results 
that originate from various values of those variational 
parameters that are not shown along the abscissa. 
With regard to panel a) it should be noted that the 
large spread of values for larger values of wh is due 
to finite accuracy: the vacuum polarization energy 
is computed by adding up several components, dis- 
played, e.g., in fig. 2] that each have a finite accuracy. 
For large values of wh the numerical uncertainty is 
of the order of the vacuum polarization energy; this 
of course contributes to the spread making it larger. 
For further comments, cf. fig|H 
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FIG. 4: Comparison of the gradient expansion results (-EDorExp, eq. 13.4311 ) with the full calculation (iJvac, eq. II3.2UI| '1. We 
also disentangle the various contributions to E^axi- The results are shown as functions of wh, and we furthermore split the 
remaining parameters into sets with ap < 0.8 (a) and ap > 0.8 (b). 




FIG. 5: The extremal (negative) binding energy, f' as a function of Nf, the number of populated bound state levels. 



In figure ini we display the binding energy eq. 14. for prescribed values of Nc and Nf, as a function of wh to 
demonstrate that we indeed obtain a local minimum in the energy when the total fermion number is constrained. 
Again, we see that a bound object requires a large value for the number of colors. 

The additional variational term ap in eq. (|2.24|) is novel in the context of strings in the standard model and it is 
worthwhile to briefly reflect on its impact on the fermion vacuum energy. The fermion zero mode exists for ^1 = § 
only for ap = 0. In addition any non-zero value of ap increases the effective fermion mass. Thus we expect largest 
binding in the small ap regime. Indeed we find that B^^f'^ is maximal for ap —> 0. It is thus a bit surprising that ap 
actually decreases the vacuum energy, cf. figure i?vac may even become negative only for sufficiently large ap. 
As a consequence, at least for objects with moderate fermion numbers, the total energy does not increase with ap 
despite that the energy eigenvalues of the individual fermion modes get larger. This effect is shown in figure 

However, for these moderate values of Nf overall binding is not yet observed. When we turn to the cases that 
exhibit binding, the system prefers ap — > because the bound states become more important as we turn up Nf. This 
effect is shown in figure 



G. Large Fermion Mass 



We have seen that stable objects can be constructed by choosing unconventional parameters. In particular, increas- 
ing the number of fermion degrees of freedom helps to bind objects with large fermion numbers, since it suppresses 
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600 h N^=7, N^20 




FIG. 6: The binding energy, eq. I|4.6^ as function of the width parameter wh for various values of Nc and Nf. 




FIG. 7: The binding energy as function of the amplitude of the charged scalar field ap for two values of the number of populated 
levels. 



the classical energy in comparison to the fermion determinant and bound state contributions. The obvious question 
is whether there are other regions in the space of model parameters that allow for stable objects. Such a search, of 
course, cannot be exhaustive and we will therefore just report on a single calculation in that direction. The most 
obvious change in parameters that enhances the role of the fermion fluctuations is an increase in the Yukawa coupling, 
/. This change increases the fermion mass and simultaneous decreases the VEV of the Higgs field in = 2 + 1 due to 
our definition eq. (|3.15|l . Since v — wd=2+i sets the scale of the classical energy, this change represents an additional 
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enhancement of the fermion vacuum contribution. 

In figure IHl we show the resulting binding energy for our choice ruf = l.STeV. As expected we do get a bound 
object for the physical value Nc = 3 and a large number of populated fermion levels. It is amusing to see that even a 



N^=3, Nj=0 



■ I ..I ..I ^^\ I' 




FIG. 8: The binding energy as function of the width of the Higgs field wh for a large fermion mass, m/ = 1.5TeV and two 
values of the number of populated levels. 

fermion number zero object (Nf = 0) acquires binding. That is, for such large a fermion mass the vacuum becomes 
unstable already for Nc = 3, at least in the one fermion loop approximation. For mf — ITOGeV we need to tune 
Nc ^ 100 to observe the same effect. Of course, we expect that there exists an intermediate value of the fermion 
mass between ITOGeV and 1.5TeV for the which this vacuum instability sets in for Nc = 3. 



V. CONCLUSIONS 



In this paper we have studied the quantum energies of static and localized W-string configurations in a = 2 + 1 
dimensional gauge theory. This is to be understood as a precursor to a full investigation of quantum energies of W/ Z- 
string configurations in the standard model in D — 'd + 1. As we know from studies on QED flux tubes, the limitation 
to D = 2 + 1 dimensions is a reliable truncation for string-type configurations when the Z) = 3 + 1 rcnormalization 
conditions are imposed and the vacuum expectation value of the Higgs field is suitably scaled. 

We have concentrated on the fermion contributions to the vacuum polarization energy as a first calculation that 
comprehensively compares that energy to the one associated with the the trivial vacuum configuration, where the 
Higgs field sits at the minimum of its potential and the gauge fields vanish. This approach is motivated by the 
observation that there exists a fermion zero-mode for field configurations that correspond to the W-string with the 
charged Higgs field being identically zero. The population of this zero mode could yield a significant energy gain. 
The incorporation of contributions due to a single fermion mode requires us to also include the full fermion vacuum 
polarization energy because it is of identical order in both the loop and large Nc expansions, where Nc counts the 
degeneracy of the fermions, e.g. the number of color degrees of freedom. We compute the effective energy as the sum 
of the classical energy, the bound state contributions and the (renormalized) fermion vacuum polarization energy. In 
the combined limit, S — > and Nc oo this approximation becomes exact at next to leading order. 

We have employed techniques that use scattering data for the computation of the vacuum polarization energy. 
These techniques have the important and outstanding feature that standard rcnormalization conditions, formulated 
in terms of Green's functions in momentum space, can be straightforwardly imposed. This is essential for a sensible 
comparison to the energy of the translationally invariant vacuum configuration and to generalize the D = 2 + 1 results 
to Z? = 3 + 1. To make efficient use of this powerful scattering theory approach, however, we have to assume that 
the fermions in the doublet are degenerate in order to obtain a partial wave expansion. Our numerical results show 
that this contribution to the total energy is small for configurations of interest, and does not destroy stability of the 
string. 

We have found energetically favored string configurations within an ansatz of variational parameters characterizing 
the Il^-string configuration. Although the true minimum likely lies outside our ansatz, its energy can only be less. 
Though we find that the fermion vacuum polarization energy reduces the effective energy, our numerical results show 
that the large classical energy carried by the pure VF-configuration cannot easily be compensated by fermion quantum 
contributions with smaller classical energies. Rather, in the search for energetically stable configurations we are led 
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to Higgs field dominated scenarios. In turn, this result reduces the importance of the zero mode. This is a significant 
result, one that runs contrary to our expectations and demonstrates that the search for stable configurations requires 
one to fully account for the fermion spectrum in the string background, rather than merely concentrating on just the 
zero mode. Furthermore, the change in character of the fields when going from the gauge field dominated to the Higgs 
field dominated scenarios is quite interesting. In the gauge field dominated case the fields have non-zero winding at 
infinity and magnetic flux. In the other case we are essentially left with only a shallow ring in coordinate space in 
which a hole is dug into the Higgs vacuum expectation value. 

In the regime of interest, where the extension of the background fields is a few times the Compton wave-length of 
the fluctuating fermion, the number of bound states essentially grows quadratically with the widths of the background 
fields. The hole in the Higgs vacuum expectation value causes the corresponding binding. Although these states are 
generally only weakly bound, they are so numerous that they can cause an object with large fermion number to 
indeed be bound. Though we did not observe such a configuration for empirically motived model parameters, we 
have seen that a doubling of the fermion degeneracy or an increase of the fermion mass provides sufficient emphasis 
on the fermion contribution to the total energy such that stable objects emerge. Typically these stable objects carry 
fermion numbers of about one hundred or even more and the configurations are quite wide, a few times the Compton 
wave-length of the free fermion. 

We have not taken the effort to determine a minimal value for the fermion mass for which a stable configuration 
emerges. However, for the empirical value of the fermion mass the critical value Nc = 7 is not too far away from 
the physical datum. It thus does not seem totally absurd to imagine that in D = 3 + 1 a stable object exists when 
adopting the standard model parameters. This certainly motivates the corresponding extension of the present study, 
but some technical obstacles must be overcome first. In 13 = 2 + 1 it was sufficient to compute Feynman diagrams 
(and Born terms) up to second order in the external fields. The expansion with respect to external fields is gauge 
variant, but only gauge invariant combinations of Higgs and vector fields have well defined Fourier transforms that 
enter the Feynman diagrams. (The full effective energy is, of course, gauge invariant as the analogous peculiarity is 
contained in the Born series for the scattering data.) Already in the present case we had to introduce a fake Higgs field 
to circumvent that problem. Once we go to higher order in the Feynman diagrams, additional manipulations of this 
kind will be needed. One possibility is to unwind the string configuration at some distant point in space by making 
the angle i^i in the parameterization, eq. H2.24|l space dependent, such that it vanishes at r ^ cxi. Unlike in the case 
of the QED flux tubes the return flux associated with 'unwinding' the strings will contribute to the energy even if 
the distant point is sent to infinity. Though that contribution is easy to estimate (and thus subtract) for the classical 
energy, it is not so for the vacuum polarization energy, which is a non-local functional of the background fields. We 
expect, however, that the artificial 'unwinding' part of the effective energy can be quantified by some approximate 
techniques such as the gradient expansion. Since the present study indicates that — > for the energetically favored 
configuration, it is also worthwhile to explore pure Higgs background configurations in £) = 3 + 1. Studies in that 
direction are in progress. 
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APPENDIX A: PERTURBATION THEORY 

1. Renormalization Conditions 

The counterterm Lagrangian for the SU{2)l model is given by 

C%'^ = citr (lU^^'lU^.) + C2 D,c^) + C3 - v^] + a [{<f>^) ~ v^]\ (Al) 

The counterterm coefficients corresponding to the renormalization conditions outlined in the main text are: 

Ci = -^/p), c2 = ^/f(eH), c:, = -fmf and C4 = ^m//f (C/^) , (A2) 
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with C4 = 4w^C4 + C3 + m^C2 and — ran/raf. For all parts of the Feynman diagram calculation we employ 
dimensional regularization. This is the reason for C3 being finite, even though it is used to cancel a linear (by power 
counting) divergence. 

The Feynman parameter integrals are given in D=2+l by: 

= (A3) 
^fHC) = l + |ln|^, (A4) 



l2\0 = 8|3(4^-(4 + ain|:p|), (A5) 
= ^(4e + (4-ain|^). (A6) 



2. The Renormalized Second Order Feynman Diagram 

The second order renormalized Feynman diagram involves only the scalar profile functions. It is most easily 
expressed using the Fourier transforms 



fnik)^ pdpMkp){l- fnip)) s.mlfp{k)= pdpJo{kp)fp{p). (A7) 
Jo Jo 

As discussed at the end of section IlII El it suffices to consider ^1 = 0. Then, the contribution to the energy is 

/>OC 

= 2m' / rjdvi Kv' + eH)mH) + 2h{iv) - 2h{£,H)\fU^}V) + W H^h) + h{iv) - AfU^i^)) (A8) 
Jo 

with — mjj /ruf and irj being a space like momentum fraction. 

3. The counterterm contribution to the energy 

In this section, we give the contribution of the counterterms to the energy that weren't used in renormalizing the 
second order Feynman diagram. The energy contribution is 

E^J'i"^ ^271 J pdpS^.t. (A9) 

with 



= -4ci sin^(^?)4# + C2^^4 - fof + fpfl) ~ c,v\{l -fl- flf 4(1 - fuf). (AlO) 
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